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ABSTRACT 

In this Letter we revisit arguments suggesting that the Bardeen-Petterson effect can coalign 
the spins of a central supermassive black hole binary accreting from a circumbinary (or cir- 
cumnuclear) gas disc. We improve on previous estimates by adding the dependence on system 
parameters, and noting that the nonlinear nature of warp propagation in a thin viscous disc af- 
fects alignment. This reduces the disc's ability to communicate the warp, and can severely 
reduce the effectiveness of disc-assisted spin alignment. We test our predictions with a Monte 
Carlo realization of random misalignments and accretion rates and we find that the outcome 
depends strongly on the spin magnitude. We estimate a generous upper limit to the proba- 
bility of alignment by making assumptions which favour it throughout. Even with these as- 
sumptions, about 40% of black holes with a > 0.5 do not have time to align with the disc. 
If the residual misalignment is not small and it is maintained down to the final coalescence 
phase this can give a powerful recoil velocity to the merged hole. Highly spinning black holes 
are thus more likely of being subject to strong recoils, the occurrence of which is currently 
debated. 
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1 INTRODUCTION 

The cosmic evolution of the spin of supermassive black holes 
(SMBH) residing in galaxy centers has several important astro- 
physical implications. By regulating the efficiency of conversion of 
matter into radiation, e, the spin magnitude a effectively determines 
the ability of black holes to grow rapidly by accretion. Defining the 
Eddington time as 
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4.5 x 10 8 yrs, 



(1) 



alescence determines the waveform of the gravitational radiation 
that may be detected by planned gravitational wave observatories 
such as NGO-eLISA. Moreover, when the spins of the two merg- 
ing black holes are significantly misaligned to each other, the rem- 
nant black hole can receiv e r ecoil velocities up to ~ 4000 km/sec 
dCampanelli et al.l 12007 jlbl lLousto & Zlochowerl 1201 ll) . There 
have been claims of the detec tion of several fast recoiling black 
holes via their emission lines jKomossaetalll 2008 ; ICivano et al.1 
20121). but the interpretation is uncle ar teogdanovic et al.| 20091 : 



Dotti et alj[2009al ; iBlecha et alj|2012h (see iKomossal |2012| for a 



where k is the Thomson opacity, c the speed of light and G the 
gravitational constant, a black hole accreting continuously at the 
Eddington rate grows in mass exponentially on a timescale r = 
efEdd/(l _ e). Rapidly spinning black holes have e as 0.5, and thus 
t ~ teddi while for non-rotating black holes r a 0.06/Edd- If most 
holes spin rapidly, the existence of luminous qu asars at redshift 
> 7, when the a ge of the Universe wa s ~ 2fEdd, dFan et al.ll2004l : 
iFan et all |20o3 : iMortlock et all 1201 ll) requires very massive (> 
10 5 M o ) seed black holes (e.g. lVolonteri & Reesll2005h. that might 
form through direct c ollapse of primordial gas jLodato & Nataraian] 
l2006h . Alternatively, iKing & Pringiel bOOd |2007|) and lKing et al.l 
1 12008!) argued instead that the spins (and so e) could remain rela- 
tively low, if most SMBH accrete through a sequence of randomly 
oriented accretion events, in what is known as the chaotic accretion 
picture. 

Black hole spin is also important in a completely different con- 
text. The relative orientation of spins during binary black hole co- 



the accretion disc 


Schnittman & Krolikl 120081: iRossi et alj|20ld: 


ICorrales etal.ll201C 


; Anderson et alj|2010h to near-Eddington lu- 



mtnosities, either in the infrared (Schnittman & Krolik 2008) or at 
high er frequencies (UV to X-rays) if the radiation is not thermal- 
ized teossi et al]20ld) . 

Black hole mergers are expected to occur in gas rich environ- 
ments, and circumbinary gas discs may play a key role in bringing 
the binary to separations close enough for gravitational wave emis- 
sion to drive the final coalescence. The effectiveness of such gas- 
assisted mergers is debated. Galaxy-scale simulations appear to 
imply that in the presence of gas a black hole pair might reach sep- 
aration s of the order of 0.01 pc and form a bina ry within 1 - 5 x 10 7 
years jEscala et alj 120051; iDotti et alJ l2009bl ). However, higher- 
resolution computations " jLodato et al.l2009h show that, at distances 
of the order of 0.01 pc, the circumbinary disc probably becomes 
self-gravitating and fragments, severely limiting the possibility of 
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mergers within a Hubble ti me. In contrast, in the c o ntext of the 
chaotic accretion picture of iKing & Pringlel d2006h . I Nixon et al.l 
d2011al) have shown that a sequence of accretion episodes where 
the circumbinary disc can be either co- or counter-aligned with 
the binary is much more effective than a simple prograde disc in 
bringing the binary to coalescence. 

Gas discs can also affect the mutual orientation of the spins 
of the two black holes through the Bardeen-Petterson (B P) effect 
jBardeen & Pettersonl Il975l ; IPapaloizou & Pringielll983l) . Lense- 
Thirring precession and viscous torques in a disc around a black 
hole tend to co- or counteralign its s pin with the disc axis 
dLodato & Pringlell2006l : Iking et alj|2005l) . In a binary black hole 
system, each hole may have an individual disc around it, whose ini- 
tial orientation agrees with that of the circumbinary disc. Then the 
BP effect can co- or counteralign the two spins. Clearly, for this to 
work, the orientation of the circumbinary disc, which is feeding the 
individual black hole discs, must stay roughl y constant. (This do es 
not hold in the chaotic accretion picture, cf iNixon et alj|20~l 1 al ibi 
lNixonll2012l.) 



Bogdanovic et al.l d2007l) (hereafter BRM) make order-of- 



magnitude estimates of the BP effect. They consider the alignment 
timescale for a single black hole and its accretion disc. BRM find 
that this is very short, and conclude that each black hole is effec- 
tively aligned with its disc and so that the two black hole spins 
are each co- or counter-aligned with the common circumnu clear 
disc at the time when they form a binary. IPerego et all 120091) im- 
prove their st udy finding s lightly longer timescale. The numerical 
simulations of lDotti et al.l J[2010h appear to support this picture, al- 
though with some differences, depending on whether the circum- 
nuclear disc is warm or cold, with the coldest (and thinnest) disc 
providing more effective alignment. 

The estimates of BRM have a number of shortcomings. First, 
as mentioned above, their argument assumes that the discs sur- 
rounding each black hole share their orientations with the circumbi- 
nary disc, and that this does not change with time. If the accretion 
process during the merger is chaotic, both assumptions may be in- 
valid. But even within the framework where the disc orientations 
track the binary orbit, the alignment timescale may not be as short 
as implied by BRM. Indeed, BRM evaluate the accretion timescale 
is by considering the flow properties at the Bondi radius, at a dis- 
tance of 40 pc from the hole. However, the viscous time at these 
distances from the hole is much larger than a Hubble time for a 
thin disc, and still larger than the binary evolution timescale if the 
disc is assumed to be thick. Thus, such estimates may not accu- 
rately describe the flow properties at the warp radius, located at a 
few hundreds S chwarzschild radii. Secondly, and most importantly, 
both BRM and lPerego eTall J2009I) assume that the warp diffusion 
coefficient, usually called y 2 , is given by v 2 /vi = I /2a 2 , where 
vi is the disc viscosity and a the standard Shakura-Sunyaev pa- 
rameter dShakura & Sunvaevlll973l). This is correct only to first or- 
der for small amplitude warps dPapaloizou" ■fePringldl 19831), while 
large amplitude warp s diffuse much less efficiently jOgilviell 19991; 



Lodato &Price]|2010|) (note also that the simulations bv lDotti et al.l 
20ld assume that the warp diffusion coefficient is independent of 



amplitude). It is thus possible that large initial misalignments stay 
misaligned throughout the merger. 

In this Letter we revisit the arguments of BRM and calculate a 
lower limit to the probability that the black hole spin remains mis- 
aligned (i.e. it does not have sufficient time to reach the aligned con- 
figuration) with its own disc at the point when gravitational wave 
reaction begins to control the merger process. In such a case it is 
very unlikely that the the two spins in a binary will end up aligned 



unless the initial misalignment is very small. We repeat that even 
in cases where we predict local alignment, a mutual coalignment of 
the two spins is only possible if the merger does not occur through 
a sequence of chaotic accretion events. 

The paper is organized as follows: in section [2] we describe 
the basic assumptions of our model and the procedure we adopt 
to compute the alignment probability. In section [3] we discuss our 
main results for a suitable choice of parameters, and in section [4] 
we draw our conclusions. 



2 MODELING SPIN ALIGNMENT 



iNataraian & Pringlel d 19981) estimate the alignment timescale for 
the BP effect in the linear case where v 2 /vi = I /2a 2 . BRM use 
these results to estimate spin alignment. In general, one can de- 
fine an a 2 coefficient corresponding to the diffusion coefficient v 2 , 
such that v 2 /vi = ail a. The value of o- 2 can be computed ba sed on 
the nonlinear theory of warp propagation of lOgilvia |l999j). This 
uses conservation laws to work out the connections between the 
viscosity coefficients, assuming that the local viscosity is isotropic. 
The validity of this th eory has been demonstrated numerically by 
lLodato& Price] J2010t) . For small warps, we have: 



a 2 _ 1 4(1 + la 1 ) 
a 



2a 2 4 + a 2 



(2) 



The a bove equation is, however, inadequate for large amplitude 
warps. lOgilviel Jl999l) provides a slowly converging Taylor series 
including additional amplitude dependent terms, but the full non- 
linear result must be computed numerically. 

The warp amplitude is measured by the parameter 
i]/ = Afil&lnR, where /? is the local inclination of the disc. To 
compute the warp evolution accurately, one would thus not only 
need to know the misalignment between the outer disc axis and the 
black hole spin, 9, but also how steep its gradient is. This clearly 
requires a detailed calculation of the shape of the disc. For our 
preliminary assessment of the effect, we simply assume that the 
warp is a smooth function, varying over a distance of the order of 
the disc size, and so approximately if/ x 9. This is a conservative 
assumption, since a steeper warp profile would increase the warp 
amplitude and thus reduce further the effectiveness of warp propa- 
gation. Indeed for sufficiently high inclination or low viscosity the 
disc m ay break, i.e. make a sharp transition between two distinct 
planes iNixon & Kingll2012iJ) . 

The alignment t imescale for small warps and a steady state 



Ihe alignment timescale tor small warps and a steady state 
disc is given by (e.g. Scheuer & Feiler|[l996l : INataraian & Pringlel 



ll998l : lLodato & Pringlel2006l) as 



3d 



v, M I R, 



V2 M \Rw 



1/2 



(3) 



where M is the black hole mass, a the hole spin parameter, M the 
accretion rate through the disc, Rs = 2GM/c 2 is the Schwarzschild 
radius, and R w is the warp radius. In <[3j the warp radius R„ is taken 
as the point at which the warp diffusion timescale equals the Lense- 
Thirring timescale, so that 



2GJ h 



(4) 



c v 2 

where J h = aGM 2 /c is the black hole angular momentum. 
The ratio R w /Rs is give n by (cf. the analogous expression in 
INataraian & Pringlel 19981) : 
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Figure 1. The value of the diffusion coefficient aj as a function of the mis- 
alignment 6 for several choices of ff=0.1 (solid line), 0.2 (dotted line), 0.3 
(dashed line), 0.4 (long-dashed line) and 0.5 (dot-dashed line). While for 
small misalignments (9 a 0), a 2 can be quite large, its value drops signifi- 
cantly (even below unity) when the misalignment approaches it/2. 



R s V' 3 UJ \R> ' 



(5) 



where H/R is the aspect ratio of the disc. Inserting this in equa- 
tion[3]and scaling the accretion rate in units of the Eddington value 
Msdii = cM/tEM we obtain: 



falign ~ 3. 37a ( — 
\a 2 



2/3 



7 x 10 6 
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\a 2 



2/3 
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M 

• lMEdd 
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(m) (crr) yr ' 



(6) 



where the various parameters have been scaled to typical values. 
Note that this timescale is only a factor 2 smaller than the typical 
estimates for th e shrinking time in a gas-rich environm ent, which 
is ~ 10 7 years tescala et alj|2005t iDotti et alj|2009bh . For small 
warps, a 2 can become very large, so the alignment timescale is 
much smaller than the shrinking time. The same also occurs for 
slowly spinning black holes, for which both a and the accretion 
efficiency e are small. Conversely, for maximally spinning black 
holes, awl, for which e » 0.4, alignment for large warps might 
requ ire a timescale longer than t he shrinking time. 

iNataraian & Pringlel d 19981) also evaluate H/R as a function 
of the main disc parameters (black hole mass and accretion rate). 
However, the dependence turns out to be rather weak, and so we 
have just chosen to keep that number fixed at a representative value 
of 0.01. In l[6j we scaled the Eddington ratio M = 0.1M EM , which 
is typical for accreting black holes in AGN. Also this assumption is 
favourable for the alignment process since it is expected that BH in 
binary systems accrete ad a reduced rate du e to binary torques (typ- 
ically at 10 per cent of the unperturbed rate. lArtvmowicz & Lubowl 
1 19961) and because frag mentation in the out er disc might lead to 
significant star formation l lLodato et alj2009l) . However, in our cal- 
culation we explore a wider range of M. 

In our calculation we take into account the relationship be- 
tween the spin parameter a and the accretion efficiency e(a), based 



on the formalism of lBardeenl d 19731) (cf. lKing et ai1l2008l , their fig. 
5): 



1 - 



2r 
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(7) 



where r isco = R isco c 2 /GM is the radius of the innermost stable orbit 
in units of the gravitational radius and is given by 



tfisco = — [3 + Z 2 + (3 - Z, ) 1/2 (3 + Z, + 2Z 2 ) 1/2 ] , (8) 
Zi = 1 + (l - a 2 )" 3 [(1 + a) 113 + (1 - a) 113 ] , (9) 

Z 2 = (3a 2 +Zf)" 2 , (10) 

and the upper/lower sign refers to particle orbiting in pro- 
grade/retrograde orbits with respect to the rotation of the hole. 

Finally, for any given misalignment 8, we compute 02 (0,6) 
numerically, using the full nonlinear theory of lOgilvia dl999h . We 
plot in Fig.[T]the resulting value of a 2 as a function of 8 for various 
choices of a, which shows that while for small misalignments (8 « 
0), a 2 can be quite large, its value drops significantly (even below 
unity) when the misalignment approaches tt/2. 



3 RESULTS 

We compute the probability of alignment through a Monte Carlo 
realization of N = 10 4 events. The initial angles 8 are randomly 
distributed between -n and n, thus including also cases where the 
disc and the black hole end up in a counteraligned configuration. 
We generate the angles uniformly in cos 8 to have an isotropic dis- 
tribution in three dimensions. For each event we generate a random 
value of the accretio n rate M/Meu betw een 10~ 4 and 1 (the same 
range as explored bv lPerego et alj2009T) . with uniform logarithmic 
distribution. For each choice of the remaining two free parameters 
a and a, we thus obtain the probability distribution of f a ij gn : we com- 
pute this distribution varying over 5 different values of a and 100 
differen t value s of a. Note that, for large 6 and small a, the theory of 
lOgilviej l ll999h would fomally predict a negative azimuthal viscos- 
ity coefficient. The behaviour of the disc in such cases is unclear, 
but since the viscosity must already have passed through small pos- 
itive values it is likely that it has in reality brok en into two dis- 
tinct planes, thus making a lignment even slower dLodato & Pried 
l20irilNixon & Kindl2012ah . When this occurs, we simply assume 
that the real alignment time becomes much larger than any other 
timescale in the problem. 

We plot in Fig.[2]with black lines the probability that the align- 
ment time is smaller than 10 Myrs (left panel) and 50 Myrs (right 
panel) as a function of the spin parameter a, for various choices 
of a (with the same notation as Fig. [T}. The values of 10 and 50 
Myrs are taken as proxies for the shrin king timescale of the binary 
jEscala et al1l2005l : IDotti et alj|2009bl) . The plots demonstrate the 
general features already discussed in section [2] For small spin pa- 
rameters alignment is indeed efficient and we expect most spins to 
align or counteralign with their discs by the time the binary ap- 
proaches coalescence. However this is not true for larger spins. In 
particular, for a > 0.5 we expect a sizeable fraction of the systems, 
of the order of 30-40% (or even more for extreme Kerr black holes), 
to not reach the aligned configuration after 10 Myrs. The result is 
still present, for 50 Myrs, in which case obviously the black holes 
have more time to become aligned. 

The probability distributions computed with a = 0. 1 are quite 
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Figure 2. Probability that the alignment time r a ij gn is smaller than 10 Myrs (left panel) or 50 Myrs (right panel), as a function of the spin parameter a, for 
various choices of a, using the same notation as in Fig.[JJ For large values of a a sizeable fraction of systems are likely to have residual spin misalignments 
after a binary shrinking time. The red lines show the corresponding probabilities in the case where we use Eq. (2) to compute t»2 and we thus do not consider 
the effect of the non-linearity of the warp on the diffusion coefficients. In this case, there is a much stronger dependence on a and for small a alignment is 
much more efficient. 



different from the others, because in this case we completely re- 
move from the alignment process all the events w ith negative vis- 
cosity coefficient. As noted above. lOgilvid JT9 991) predicts this be- 
havior for small values of a. We note however that even such a 
strong assumption does not affect significantly the distributions for 
large value of a where the effects described above are definitely 
more important. 

For all the other values of a, there is no negative viscosity 
event: we note that our results are almost independent on a and 
that they are not very sensitive to the proxy chosen (10 or 50 Myrs), 
ensuring that our results are robust with respect to different choices 
of the merger timescale. 

In Fig. [2] we also plot with red lines (and using the same line 
styles as above) the corresponding probabilities in the case where 
we use the warp diffusion coefficien t in the linear regime (Eq. 2) 
dQgilvid[l999l ; lLodato & Pricejl2010h and we thus do not consider 
the effect of the non-linearity of the warp on the diffusion coeffi- 
cients. The effect of the non-linearity is clearly very significant. In 
particular, it turns out that in the linear case the alignment proba- 
bility has a much stronger dependence on a (essentially because in 
the non-linear case the range of values of or 2 is much smaller than 
in the linear case) and alignment is indeed much faster, especially 
for low a. 

We have also considered simpler cases where we keep the Ed- 
dington ratio fixed, in order to evaluate more clearly the depen- 
dance on this parameter. Clearly, since the alignment time is in- 
versely proportional to M, larger Eddington ratios (such as those 
assumed by BRM) imply a faster alignment. However, also in 
these cases, we find that alignment within lOMyrs occurs only for 
a < 0.6 and a < 0.3 for the two cases where M = 0. IMecM and 
M = 0.01M Edd , respectively. In order to have fast alignment for 
high spins the accretion rate is required to stay at the Eddington 
level for a prolonged period of time. 



4 DISCUSSION AND CONCLUSIONS 

We have revisited the arguments suggesting that gas disc are ef- 
fective in bringing the spins of the two black holes in a merging 
binary into alignment. In particular, we have improved on previous 
estimates (BRM) by taking into full consideration how the align- 
ment timescale changes with the system parameters. In addition, 
we have also included the reduction in the warp diffusion coeffi - 
cient when the misalignment angle becomes large togilvid[l999h . 
an effect previously neglected. 

Contrary to previous claims, if the black holes are rapidly 
spinning (a > 0.5), the system is would not end up to be com- 
pletely aligned, in up to 40% of the cases, at the time at which the 
two holes are brought together at distances of the order of 0.01 pc 
(the current resolution limit for numerical simulations of the pro- 
cess). 

Our results are consistent, at l east within orders o f magnitude, 
with the previous investigation bv lPerego et al. I d2009h . In particu- 
lar, we find a similar dependence of the alignment timescale with 
respect to the spin of the bl ack hole and the Eddington ratio. Con- 
trary to lPerego et al.l d2009l) , our estimates show that the alignment 
timescale is a strong function of the initial misalignment angle. This 
comes from having considered the full non linear warp propagation 
theory, instead of restricting to the small amplitude regime: this 
is a key element, because it allow us to compute a probability for 
the alignment process, based on the expected distribution of initial 
misalignments. 

Highly spinning black holes are thus more likely to maintain a 
significant misalignment, which causes a high recoil velocity. The 
current lack of observational evidence for strongly recoiling black 
holes (but see the recent case of CID-42, Icivano et alj|2012l) sug- 
gests that the average black hole spin is rat her low, in line with th e 
predictions of the chaotic accretion picture <King & PringiefeoOfj) . 
We note that for a < 0.6, as suggested above, the accretion effi- 
ciency is e =! 0.1. This is in line with the values sugg ested by the 
Soltan argument JSortarJll982t l Yu & Tremaindl200 Z) relating the 
average SMBH mass to the radiation background of the universe 
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(see also detailed discussion in iKing et alj[2008h - We stress here 
that our conclusion regarding the magnitude of the black hole spin 
is related to the average properties of the BH population. Individ- 
ual SMBH mig ht well have large spins, and o bservations of broad 
iron lines (e.g.. lBrenneman & Revnoldsll2006h would naturally be 
biased in their favour. On the contrary, if recoiling black holes are 
found to be more common, such limitation on the magnitude of the 
black hole spin would not apply. 

Our model is certainly very idealized and could be refined in 
several ways. First of all, we have used a uniform logarithmic distri- 
bution of the Eddington value withou t considering in a detailed way 
how do tidal effects and gap opening dArtvmowicz & Lubowl 19961) 
affect the accretion rate on the binary elements. Another interest- 
ing effect arises from the fact that the alignment timescale scales 
inversely with the accretion rate. In binaries with a large mass ra- 
tio, where accretion occurs preferentially onto the secondary, the 
primary might be much harder to align, with potentially interesting 
effects on gravitational wave emission. 

A second important limitation of our work comes from the 
assumption that the lengthscale over which the disc inclination 
varies is comparable with the disc size. Also in this case, our as- 
sumptions go in the direction of underestimating the alignment 
timescale, since if the warp occurs over a short lengthscale even 
a relatively small initial misalignment might be more difficult to 
realign. To take this effect into account, one would need to ex- 
plici tly solve for the disc sh ape in a time-dependent calculation 
(e.g. lLodato&Pringlell2006l) . As a consequence, we cannot quan- 
tify the exact degree of misalignment at the end of the merger, but 
can only perform a simple comparison between the two timescales. 
In order to quantify more precisely the level of misalignment, a de- 
tailed ti me-dependen t calcu latio n should be perfor med, along the 
lines of IPerego et"aTI J20091) and lDotti et all feOld) . Additionally, 
our choice of 10 and 50 Myrs as proxies for the shrinking timescale 
is clearly very simplified, and might in particular be a function of 
the system parameters, such as the accretion rate. For example, a 
large value of f a i ign can be due to a small value of M and should 
then be compared with a longer shrinking timescale, suitable for 
the expected gas-poor environment. 

Third, we have considered only the effect of gas discs on the 
evolution of black hole spin. To follow the evolution of the spins 
after decoupling due to relativistic effects, a post-Newtonian ap- 
proach is ne eded, taking into account the possib le role of spin-orbit 
resonances dSchnittm an 2004; Berti et alj|20121) . 

Finally, we have considered only the alignment process of one 
black hole with its own disc. Aligment here is clearly a necessary 
condition for coalignment of the two spins in a merging binary, 
Even in cases where we predict alignment of individual spins and 
discs, the coalignment of the two spins with the larger scale cir- 
cumbinary disc can present add itional severe difficultie s if the black 
holes are rapidly spinning (see iNixon & Kingll2012bl for a discus- 
sion). All these effects act in the direction of reducing the likelihood 
of alignment of both spins even further. 
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